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Abstract 

We study the time evolution of a periodically driven quantum- 
mechanical system coupled to several reservoirs of free fermions at 
different temperatures. This is a paradigm of a cyclic thermodynamic 
process. We introduce the notion of a Floquet Liouvillean as the gen- 
erator of the dynamics of the coupled system on an extended Hilbert 
space. We show that the time-periodic state which the state of the 
coupled system converges to after very many periods corresponds to a 
zero-energy resonance of the Floquet Liouvillean. We then show that 
the entropy production per cycle is (strictly) positive, a property that 
implies Carnot's formulation of the second law of thermodynamics. 



1 Introduction 



During the past several years, there has been substantial progress in the 
program of deriving the fundamental laws of thermodynamics from nonequi- 
librium quantum statistical mechanics (see [A-SFl, A-S] for a synopsis). In 
this paper, we make a contribution to this program by studying Carnot's 
formulation of the second law of thermodynamics from the point of view of 
quantum statistical mechanics. For the sake of concreteness, we consider 
a periodically driven two-level quantum- mechanical system, S, coupled to 

* e-mail: walid@itp.phys.ethz.ch, juerg@itp.phys.ethz.ch 
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n > 2 reservoirs of free fermions, TZi, ■ ■ ■ ,TZn- Our analysis can be general- 
ized to a system composed of an arbitrary quantum-mechanical system with a 
finite-dimensional Hilbert space coupled to several reservoirs of (free) bosons 
or fermions at different temperatures.^ In order to study the time evolu- 
tion of the coupled system, we extend Floquet theory for periodically driven 
quantum systems at zero temperature (see for example [Ho, Yal, Ya2]) to 
apply to systems at positive temperatures. 

In particular, we introduce an operator, the Floquet Liouvillean, gener- 
ating the dynamics of the system on an extended Hilbert space, and we show 
that the time-periodic state, which the state of the system converges to after 
very many periods, corresponds to a zero-energy resonance of the Floquet 
Liouvillean. We also prove (strict) positivity of entropy production per cycle, 
which amounts to Carnot's formulation of the second law of thermodynamics. 
For weak enough coupling of the small system, S, to the reservoirs, the time- 
periodic state corresponding to the zero-energy resonance can be expanded 
in powers of the coupling constant. This is of considerable practical impor- 
tance, since it enables one to explicitly compute the degree of efficiency of 
the coupled system. Further discussion of the second law of thermodynamics 
and another proof of convergence to time-periodic states using methods of 
scattering theory will appear in [A-SFl]. As far as we know, all previous 
investigations of Carnot's formulation of the second law of thermodynamics 
from the point of view of quantum statistical mechanics assumed that the 
coupling is switched on at time t = to, and then switched off at a later time 
t = to + T (see [PW, BR] and references therein). The novelty of our approach 
is to prove that the state of the coupled system converges to a time-periodic 
state with the same period as the one of the interaction, and that entropy 
production is (strictly) positive; (see also [FMSUe,MMSl,2]). 

The organization of this paper is as follows. In section 2, we recall 
some basic notions from quantum statistical mechanics, in particular, time- 
dependent perturbations of C*-dynamical systems, relative entropy, and Carnot's 
formulation of the second law of thermodynamics. We also discuss sufficient 
conditions to prove strict positivity of entropy production per cycle. These 
conditions are satisfied in the concrete example we consider in the following 
section. In section 3, we discuss a concrete model, and we state the as- 

"'^Bosonic reservoirs are more difBcult, technically, since the interaction term coupling 
E to the bosonic reservoirs is generally an unbounded operator. However, one may readily 
extend the methods developed in [MMS1,2] to study the bosonic case. 



2 



sumptions we make on the interaction between the small system, S, and the 
reservoirs, TZi,--- ,7ln- In section 4, we introduce the Floquet Liouvillean, 
whose spectrum we study using complex spectral deformation techniques. In 
section 5, we use results on the Floquet Liouvillean to prove convergence of 
the state of the coupled system to a time-periodic state. This is one of the 
main results of our paper. We then prove strict positivity of entropy pro- 
duction per cycle in section 6, and we discuss how to compute the degree of 
efficiency of the coupled system for weak enough coupling. The main ingre- 
dients of our analysis are a concrete representation of the fermionic reservoirs 
(Araki-Wyss representation), a spectral approach to cyclic thermodynamic 
processes using the so called Floquet Liouvillean, and complex spectral de- 
formation techniques. 

Acknowledgements W.A.S. would like to thank Gian Michele Graf for 
enjoyable discussions. 

2 General considerations 

In this section, we recall some of the basic notions of quantum statistical 
mechanics, and we introduce the notion of time-periodic states. Although 
some of the material is standard, it is presented in this section in order to 
make our exposition reasonably self-contained. 

In the algebraic formulation of quantum statistical mechanics, a physi- 
cal system is described by a C* or W*— dynamical system. Since we only 
consider fermionic reservoirs in this paper, we restrict our attention to the 
discussion of C*— dynamical systems. However, our analysis can be general- 
ized to VT*— dynamical systems; (see for example [DJP], and also [MMS1,2]). 

A C*-dynamical system is a pair (0, a), where is the kinematical alge- 
bra of the system, a C*-algebra with identity, and a, which specifies the time 
evolution, is a norm-continuous one-parameter group of ^-automorphisms of 
0. A physical state of the system is described by a positive, linear functional 
uj, with a;(l) = 1. The set E{0) of all states is a convex, weak-* compact 
subset of the dual 0*. 

Physically relevant states of (isolated) thermal reservoirs are assumed to 
be normal to equilibrium states characterized by the Kubo-Martin-Schwinger 
(KMS) condition. An equilibrium state of (0, a) at inverse temperature 
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/3, u;^, is an (a,/9)-KMS state satisfying 

for all a, 6 G 0", where 0° is norm-dense in 0. 

We briefly recall the perturbation theory for C*— dynamical systems; (for 
further details, see for example [BR]). Let S be the generator of a, ie, a* = 
e**^, t G R. The domain of the derivation 6, ^{6), is a *-subalgebra of 0, and 
for all a,beV{5), 

S{ay = 6{a*), 6{ab) = 6{a)b + a6{b) . 

Consider the time-dependent family of perturbations, {gV{t)}te'R, with 
selfadjoint elements gV{t) G 0. Then ag, the perturbed time evolution, is a 
norm-continuous one-parameter family of ^-automorphisms of satisfying 

f^al{a)=al{5ia) + tg[Vit),a]), (1) 

and a°(a) = a, for all a G 0. Explicitly, 



a*(a) = a*(a 



/•I Ptl Ptn-l 

(2) 



n>l 



In the standard interaction picture. 



alia) = Tla\a)Tl*, 



where Tg is a unitary element of which satisfies 



d 
Jt 



-Jl = K^\9Vm (3) 



and = 1 ; ie. 
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n>l 



Next, we discuss the notion of relative entropy and of entropy produc- 
tion. Assume that there exists a reference C*-dynamics on and a state 
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uo with the property that uo is an (cr^j, — 1)-KMS state. (Equivalently, at in- 
verse temperature /5 > 0,cr^^ = aj*^^). Let 8^ be the generator of cr^, and 
let {Tiui, TTt^, f^oj) be the GNS representation of the kinematical algebra as- 
sociated to the state uj. (For further discussion of the GNS construction see, 
for example, [BR].) 

A state 7] G E{@) is called (^-normal if there exists a density matrix 
on 7i(^, such that, for all a G 0, 

7]{a) = Tr^prjTTujia)), 

where Tr is the trace over Ticj. We will denote by A/L the set of all cu-normal 
states in E{0). 

For a state rj G A^,, which might be time- dependent, denote by Ent{ri\u) 
the relative entropy of Araki [Ar,D]. For finite systems, 

Ent{ri\uj) = —Tr{r]loguj — 77 log?]). (4) 

If ?7 ^ AC, set Ent{r]\Lj) = +00.^ 

For a self-adoint perturbation gV(t) G ^{6^) of the dynamical system, as 
discussed above, we define the rate of entropy production in a state 77 G E{0) 
relative to a reference state u as 

Epiv):=v{5U9Vm: (5) 

see for example [BR,JP4]. 

It is instructive to see how one obtains this expression for entropy produc- 
tion as the thermodynamic limit of quantities referring to finitely extended 
reservoirs. [A-SF2] Consider a quantum system composed of a small system, 
S, coupled to n reservoirs, T^i, ■ ■ ■ ,7ln, at inverse temperatures Pi, - ■ ■ , Pn, 
respectively. We first treat the reservoirs as finitely extended systems and 
then take the thermodynamic limit of suitable quantities. The total Hamil- 
tonian of the finite coupled system is 

n 

H{t) = H^ + Y,H'^'+9V{t), 

1=1 

where is the Hamiltonian of the uncoupled small system, is the 
Hamiltonian of the i^^ uncoupled reservoir, and gV{t) is the interaction term 

^Note the choice of the sign of relative entropy. 
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coupling S to the reservoirs. Let the reference state of the reservoirs be uP^ 
with corresponding density matrix given by 



where is the density matrix corresponding to the equihbrium state of the 
reservoir TZi at inverse temperature which is given by 



where Tr-ji^ is the trace over the Hilbert space of the reservoir 7?.j. 

Define cuj := o a* , where uj is the initial [a^^, — 1)-KMS state, and let 
be the density matrix corresponding to lo. For a finite system, 

Ent{yjt\iJ^) = -Tr(p'^' logp^) + Tr(p'"' logp"^*) 
= -Tr(p'"' log p^) +Tr(p'^ logp^). 



and hence 



^Entiujtluj'^) = 2Tr([i7(t),p-']logp^) 



dt 

By cyclicity of the trace we have that 



7 n 

1=1 



I 

i=l 

= uoal{5^{gV{t))). 

where we have used in the last equation that 5^^ = —Yl^=iPi^iy 6i = 
i[H'^\ ■]. Note that the thermodynamic limit of the entropy production rate 
is well-defined. 

One may relate the entropy production to the heat flux from the reser- 
voirs. The heat flux from reservoir 7^,;, z = 1, ■ ■ ■ , n, at time t is 



:= -~^uoal[H-^^) (6) 

= -zu{al{\gV{t),H^^])) 

= uoa\{6,{gVm- (7) 
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It follows that 

n 

J2f^i^'^'it) = -Ep{^°(^l), (8) 

i=l 

which, for n = 1 and reversible processes, is a familiar equation. 

We now give a definition of a time-periodic state in case of a time-periodic 
coupling. 

Time-periodic state. ^ Assume that the perturbation gV{t) is time periodic 
with period r, and is norm-dijferentiable, for t > 0. For s G [0, r), define the 
time-periodic state as 



hm uj o aj^'. (9) 



" n— >oo ^ 

Note that it follows from this definition of a time-periodic state that. 

We will show in section 5 that this state is related to a zero-energy resonance 
of the so called Floquet Liouvillean. 

Next, we exhibit a connection to Carnot's formulation of the second law 
of thermodynamics. Consider a cyclic thermodynamic process in which S is 
coupled to two reservoirs, IZi and 7^2, at temperatures Ti and T2, respectively, 
with Ti > T2. Reservoir 1 acts as a heat source and reservoir 2 as a heat sink. 
Recall that the generator of the free dynamics of reservoir TZi is 5i,i = 1,2. 
Since V{t) is norm differentiable, for t > 0, it follows that 



for i = 1,2. Using the fact that z = 1, 2, is a ^-derivation which commutes 
with a* and equation Q, it follows that 

j^uj{Tl6,,{Tl*))=zujoal{6,{gVm- 



weaker definition of a time-periodic state is 
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LU^ ,. := lim — / dtto o 
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Therefore, the heat energy flowing from reservoir TZi into system S during 
the time interval [0, t] is 

AlQ, = -j'dt'uoa'^{6,{gV{t'))) (10) 

= -2^(r*5,(r**)) , (11) 

i = 1,2. Moreover, integrating (jS)), we have 

PiAlQi + = -Ent{uj o al\uj) < 0, (12) 

since the relative entropy Ent{uj o a* |u;) > 0.^ 

Now deflne the heat flow per cycle from each reservoir TZi into E as 

AQi ■= lim [Q,((n + l)r) - Qi{nT)ii = 1, 2. (13) 



n— >oo 



We assume that the system converges to a time-periodic state, and that 
during every cycle, it performs work, ie, 

AA = AQi + AQ2 > 0. (14) 

It follows from the deflnite sign of relative entropy and the existence of the 
time-periodic limit that the entropy production per cycle is nonnegative, 

AEnt = [ dtuj+,{6MVm = -iPiAQ, + P2AQ2) > 0. (15) 
Jo 

The fact that (3i < /?2, dH)) and (USD, imply that AQi > 0. It then follows 
that the degree of efficiency, 

AA , , 

AQi + AQ2 ^^^^ 



AQ 



< — =: r/^'^™"*, (18) 

Ti 

which is nothing but Carnot's formulation of the second law of thermody- 
namics. In certain situations, one can show that 

lim [Ent{uj o a("+i)^|cj) - Ent{LU o a^jcj)] > 0, (19) 



''^This follows from a general trace inequality, see for example [BR] , 
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which holds for the model we consider in this paper, and hence rj < rp^-"^™^ _ 
The following proposition states sufficient conditions for inequality (|19p to 
hold. 



Proposition 2.1 

For t G R^, let s := t mod t. Suppose that 

(a) cu^g ^ A/L, and that 

(h) sup^,R+ I dt{u:l,^,,^{5^{gVm - u; o a\{5^{gVm}\ < C, where 
C is a finite, nonnegative constant. 

Then Ep{uj+^) > 0. 

Proof. Suppose that Ep{ijjg J) = 0. Then 

Ent{uj o al\uj) = [ dt'u o a'g{6^{gV{t'))) 
Jo 

= fdt'iu o a';{6MV{t'))) - u;^,,r,,,AU9V{t'm 
Jo 

< C. 

In particular, 

Ent{u+\u) = lim Ent{uj o < C. 

^' n— >oo ^ 

Let A4 = 7r^(0)", the double commutant of tTuj{&), and let A^* be its predual. 
The set of all states 7 G A/L such that Ent{'y\uj) < C is a ( A/i *, -compact 
(see [BR,D]). It follows that G A/L, which contradicts assumption (a). □ 

Summary of main results 

Before specifying the concrete model we study, we briefly describe the 
main results of this paper, deferring precise statements and proofs to sub- 
sequent sections. One key result of this paper is transposing the problem 
of proving convergence to a time-periodic state to a spectral problem by 
introducing the so called Floquet Liouvillean (Section 4). In Section 5, The- 
orem 5.1, we show that the time-periodic state to which the state of the 
coupled system converges after very many periods is related to a zero-energy 
resonance of the Floquet Liouvillean. We also establish strict positivity of 
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entropy production per cycle in the time-periodic state (Section 6, Theorem 
6.3). In the case of two reservoirs, positivity of entropy production imphes 
Carnot 's formulation of the second law of thermodynamics, 



Our analysis also has some quantitative implications: At weak coupling, 
the time-periodic state is analytic in the coupling constant, and hence one 
can calculate the entropy production per cycle perturbatively. This leads 
to a perturbative calculation of the degree of efficiency when the system is 
operated heat engine. 



3 The Model 

As an example, we consider a two-level quantum system S coupled to n 
reservoirs ,7^i,--- ,Tln,n > 2, of free fermions in thermal equilibrium at 
inverse temperatures Pi, - ■ ■ , and chemical potentials /^i, • • • , /in- ^ 



3.1 The small system 

The kinematical algebra of S is 0^ = Ai^C"^), the algebra of complex 2x2 
matrices over the Hilbert space = C^. Its Hamiltonian is given by 
= uJoas, where ai,i = 1,2,3, are the Pauli matrices. When the system 
S is not coupled to the reservoirs, its dynamics in the Heisenberg picture is 
given by 

C(^[aj := e ae , [ZO) 

for a G 0^. 

A physical state of the small system is described by a density matrix pg. 

1 /2 

The operator = belongs to the space of Hilbert-Schmidt operators, 
which is isomorphic to ® 7i^. Two commuting representations of 0^ on 
® are given by 

7rs(a) := a® 1^ , (21) 
7r*{a) := 1^ ® C^aC^ , (22) 

^For the sake of simplicity of exposition, we set the chemical potentials of the reservoirs 
to be equal in the subsequent sections. 
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where is an antiunitary involution on corresponding to complex con- 
jugation; (see for example [BFS]). 

The generator of the free dynamics on the Hilbert space T-F ® T-F is the 
standard Liouvillean 

= -1^ . (23) 

The spectrum of is cr(L^) = {— 2ci;o5 0, 2u;o}, with double degeneracy 
at zero. 

Let uj^ be the initial state of the small system E, with corresponding 
vector Vf' G ® 7i^. The modular operator associated with is = 
up ® uj"^ , and the modular conjugation operator, J^, is given by 

for (p,ip E Tip. If tus corresponds to the trace state, then = 1^ ® 1^. 
3.2 The reservoirs 

Each thermal reservoir is formed of free fermions. It is infinitely extended 
and dispersive. We assume that the Hilbert space of a single fermion is 
h = L^(R+, m(M)(iM; where B is an auxiliary Hilbert space, and m{u)du 
is a measure on R"*". We also assume that the single-fermion Hamiltonian, 
h., corresponds to the operator of multiplication by n G R^. For instance, 
for reservoirs formed of nonrelativistic fermions in R^, the auxiliary Hilbert 
space E is L'^{S'^,da), where S'^ is the unit sphere in li^, da is the uniform 
measure on S^, and u = |/cp, where k G R^ is the particle's momentum. In 
the latter case, the measure on R"*" is choosen to be m{u)du = \^/udu. For 
the sake of concreteness, we will consider B = L"^ {S"^"^ , da) , d > 2, in the 
sequel. 

Let b and b* be the annihilation-and creation operators on the Fermionic 
Fock space J^{L'^{R+,B)). They satisfy the CAR 

{b*{f),b*{g)} = 0, (24) 
{b{f),b*{g)} = {f,g)l, (25) 

where 6* stands for b or b*, f,g E L'^ ; B) , and (■, ■) denotes the scalar 
product in L^(R"'"; B). Moreover, let fl^ denote the vacuum state in JF(L^(R"'"; B)). 
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The kinematical algebra, 0^% of the i^^ reservoir 7li,i = I,-- - ,n, is 
generated by bf and the identity The free dynamics of each reservoir 
(before the systems are coupled) is given by 

c^nMif))=bf{e''V), (26) 

for i = I,-- - ,n,f G L^(Ii~^;B). For a nonzero chemical potential, /ij, of 
reservoir TZi, an auxiliary free dynamics is generated by H^^ = dVi{h — fJ^i); 
see for example [BR]. 

The (a^,, /ij)-KMS state, of each reservoir 7li,i = I,-- - ,n, at 
inverse temperature Pi and chemical potential fit, is the gauge invariant, 
quasi-free state uniquely determined by the two-point function 

where Pft ,;.,(«) := 

Next, we introduce J^f^ := T'^^{L\R+; B)) ^ J^'^' {L\R+ ; B)), the GNS 
Hilbert space for the Araki-Wyss representation of each fermionic reser- 
voir TZi associated with the state cu^', [ArWy]. Let hi and b* denote the 
annihilation- and creation operators on J-'^' {L"^ (IV' ; B)) satisfying the CAR, 
and denote by fi^' the vacuum state in JF'^'(L^(R"'"; i3)), with fejfi^' = 0. The 
Araki-Wyss representation, vTj, of the kinematical algebra 0^% i = 1, ■ ■ ■ , n, 
on J^f^^ is given by 

7^^{h{f)) := hi^l - pft,^, /) ® + (-1)^^ ® m^yp^J) , (28) 

where Ni = (irj(l) is the particle number operator for reservoir TZi. Further- 
more, fi^' ® f^^' G JF/^^ corresponds to the equilibrium KMS state u'^* of 
reservoir TZi. 

The free dynamics on the GNS Hilbert space J-'f^ of each reservoir TZi is 
generated by the standard Liouvillean The modular operator associated 
with {0^\uj^'^) is given by 

and the modular conjugation is given by 

J^«(^ ® $) = (_l)^>W-l)/2$ ^ (_l)^.W-l)/2^^ 
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for ^, $ G J^f^] (see, for example, [BR]). 

In order to apply the complex translation method developed in [JP1,2,3], 
we map Tf^ = J^^»(L2(R+; B)) ® J^'^'{L\R+; B)) to ^^'(^^(R; B)) using 
the isomorphism between L^(R+;;B) © L^(R+;S) and L'^(R]B). To every 
/ G L^(R"''; B), we associate functions fjs^^, G Iv^(R; iS) by setting 

/,,Ka) ^^^^/ggVl -P/3>) /K^)^ ->0 ^ ^29) 
I A/m(-n)^p^,^(-u) /(-n,o-) , n<0 



and 



^''^ ' [ v^m(-n)zA/l - P/3,^(-?i) 7(-n,a),u < 

= ^7/3,m(-«'^)' (30) 

where m{u)du is the measure on R"*", see eq. PHj) . (For a discussion of this 
map, see Appendix A.) 

Let and a* be the annihilation and creation operators on JF^» (L^(R, du\ B)). 
Then 

vr.(&f(/))-af(/ft,,J, (31) 

vrf(&f(/))-^(-l)^'af(/|,,,); (32) 
^f]-^. ^ (33) 

where af stands for aj or a*, and Vt'^^ is the vacuum state in JF^*(L^(R, S)). 
Using eqs. (j77|) and one readily verifies that 

(^''%«*(/ft,MJ«.(/ft,MJ^^o = (^^%«r(/L>.(/L.)^^o 

= ..^»(6*(/)6,(/)) 

= (/,Pft,M,(-)/) • 

Moreover, the free Liouvillean on JF''^'(L^(R; i3)) for the reservoir T^j is 
mapped to 

L^' = dT,{u,) , (34) 

where G R. 
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3.3 The coupled system 

The kinematical algebra of the total system, S V TZi V • • ■ V TZn, is given by 

= 02 ® 0^1 ® . . . ® 0^" , (35) 

and the Heisenberg-picture dynamics of the uncoupled system is given by 

"o = "E ® ® • ■ ■ ® «k • (36) 

The representation oi on H := ^ T^^{L^{R; B)) ■ ■ ■ ® 
:F'^"(L^(R; B)), determined by the initial state 

uj = uj^ ^00^' u'^" (37) 

by the GNS construction, is given by 

TT = TTe ® TT^, O ■ ■ ■ ® 7r^„, (38) 
and an anti-representation commuting with tt by 

7r# = 4 ® vr* ® ■ ■ ■ ® 7r£ . (39) 

Moreover, let Vl := Q^^Cl^^- ■ ■(^Cl'^" denote the vector in Ti. corresponding 
to the state u. Denote the double commutant of 7r(0) by M. := tt{0)", which 
is the smallest von Neumann algebra containing 7r(0). 

The Liouvillean of the total uncoupled system is given by 

n 

Lo = + ^ L^^ . (40) 

i=l 

This defines a selfadjoint operator on Ti. 

For a G 0, we abreviate 7r(a) by a whenever there is no danger of confu- 
sion. The modular operator of the total system is 

A = ® A^i ® ■ ■ ■ ® A^", 

and the modular conjugation is 
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According to Tomita-Takesaki theory, 

JMJ = M\ A^'MA^'' = M , 

for t G R; (see for example [BR]). Furthermore, for a & Ai, 

JA^/'^aQ = a*Q. (41) 

The system S is coupled to the reservoirs TZi, ■ ■ ■ , 7^„, through an inter- 
action gV(t), where V{t) G is given by 

n 

V{t) = ® bUMt)) + a+ ® kim)} , (42) 

i=l 

where a± = ai ± icr2, and fi(t) G L^(R+; B),i = 1, ■ ■ ■ ,n, are form factors.^ 
The standard Liouvillean of the interacting system acting on the GNS 
Hilbert space Ti is given by 

Lg{t) = Lo + gI{t), (43) 

where the unperturbed Liouvillean is defined in ()40|) . and the interaction 
Liouvillean determined by the operator V{t) is given by 

I{t) = {V{t) - JV{t)J} 

n 
i=l 

- zl^ ® a_ ® (-l)^^a*(/4,,^(t)) - zl^ ® a+ ® i-ir^a^if^^Jt))}, (44) 

where ai,a* are the annihilation and creation operators on the fermionic 
Fock space J-'^^ {L"^ (R; B)). Note that since the perturbation is bounded, the 
domain of Lg{t) is ViLg{t)) = V{Lo). 

Let Ug be the propagator generated by the standard Liouvillean. It sat- 
isfies 

ddIg{t) = -^^git)Ugit);Ug{0) = l. (45) 

The Heisenberg-picture evolution is given by 

= U;(^)aUg{t) , (46) 

^Note that this form of an interaction preserves the total number of fermions in S V 
TZi V • • • V TZni as required by gauge invariance of the first kind. 
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for a e 0. 

Generally, the kernel of Lg(t),Ker Lg, is expected to be empty when at 
least two of the reservoirs have different temperatures; (see section 6 and 
[JP3,MMS1,2])^. This motivates introducing the so called C-Liouvillean, Lg 
, which generates dynamics on a Banach space contained in Ti (isomorphic 
to 0) and which, by construction, has a non-trivial kernel. 

Consider the Banach space 

C{0,Q) := {aQ:ae 0}, 

with norm ||afi||oo = Since Q is separating for 0, the norm ||ar2||oo is 
well-defined, and since fl is cychc for 0, C(0, fl) is dense in H. 
We set 0(a) = afl, and define a propagator Ug(t,t') by 

<f){af{a)) = Ug{t,t')an. (47) 

Then 

dtUg{t, t') = iLg{t)Ug{t, t' ) WUk Ug{t, t) = 1, (48) 

and 

Ug{t',t)n = n. (49) 

Differentiating (jTTj) with respect to t, setting t = t', and using PHj) . 
and (jHj), one obtains 



[{L, + gV{t))a-a{Lo + gVmn = 



Hence, the C-Liouvillean is given by 



[{Lo + gV{t))a-{V{t)a*y]n 

(Lo + gVit) - ^JAi/V(t)A-i/V)afi 

Lg{t)aQ . 



Lg{t) := Lo + gV{t) - gJA'/'V{t)A-'/'j . (50) 
Note that, by construction, 

Lg{t)VL = 0, 

for all t e R. 

^This is consistent with the fact that the coupled system is not expected to possess the 
property of return to equiUbrium if the reservoirs have different temperatures (or chemical 
potentials). One can verify that, indeed, this is the case when assumptions (Af)-(A3), 
below, are satisfied. 
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Next, we discuss the assumptions on the interaction. For 5 > 0, we define 
the strips in the complex plane 

■={zeC: \Imz\ < 6} 

and 

r{5) := {z e C : -5 < Imz < 0}. (51) 
Moreover, for every function / G L'^(R~^; we define a function / by setting 



f{Ha),u<0 ' ^ ^ 

where m{u)du is the measure on R+. Denote by H^{6,B) the Hardy class of 
analytic functions 

h : I{S) B, 

with 

\\h\\H\s,B) ■= sup / \\h{u + i9)\\ldu < oo. 

|9|«5 JR 

We require the following basic assumptions on the interaction term. 



(Al) Periodicity. 

The interaction term V{t) is periodic with (a minimal) period r < oo: 

V{t) = V{t + T). 

(A2) Regularity of the form factors. 

Assume that 36 > 0, independent of t and i G {1, ■ ■ ■ , n}, such that 

e-"^"/Vft,..Kt)ei/^(5,S), 

where fi is as in eq. (jUJ and fp^^ is defined in 

(A3) Fermi Golden Rule. 
We assume that 

n 

J]||/,(2a;o,t)|||>0, 

i=l 
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where fi is defined in (I52j) . which means that the couphng of S to the 
reservoirs is non-vanishing in second order perturbation theory.^ 

Let / be the Fourier transform of / given by 

Uu, uj) = - r dte-'^^^'fiu, t), (53) 
^ Jo 

where u = and r is the period of the interaction term (see assumption 
(Al)). Then f{u,t) = Emez I* follows from (A2) and Parse- 

val's theorem that 

^ \\fm{u + mLJ,Lj)\\l < oo, (54) 
meZ 

for M G R. 

Let Ug be the propagator generated by the adjoint of the C-Liouvillean, 

ie, 

dtUg{t,t') = -iUg{t,t')L;{t), (55) 
Ug{t,t) = 1. (56) 

Assumption (A2) implies that the perturbation is bounded, and hence the 
domain of Lf, where Lf stands for Lg or L*, is 

P(L#)=P(Lo), 

and Ug, Ug are bounded and strongly continuous in t and t' . 

4 The Floquet Liouvillean 

In this section, we extend Floquet theory for periodically driven quantum 
systems at zero temperature to a theory for systems at positive temperatures. 
The goal is to investigate whether the state of the coupled system converges to 

®For instance, when the reservoirs arc formed of nonrclativistic fcrmions in R'^, an 
example of a form factor satisfying assumptions (A1)-(A3) is given by 

U{u,t) = h,{t)\uY'e-\-^\\ 
where hi{t) is a bounded, periodic function of t G R and n > 2. 
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a time-periodic state. We introduce the Floquet Liouvillean, which generates 
the dynamics on a suitable Banach space, and we show in the following 
section that the time-periodic state to which the state of the system converges 
after very many periods is related to a zero-energy resonance of the Floquet 
Liouvillean. 

We consider the extended Hilbert space li, := L^([0,r]) ® H, where r is 
the period of the perturbation appearing in (j42p and (Al), and we introduce 
the Floquet Liouvillean 

K;:=-idt®l + l^L*g{t), (57) 

with periodic boundary conditions in t. Note that, under assumption (A2), 
K* is a closed operator with domain V = V{idt ® 1) fl V{1 ® Lq). 
By Fourier transformation, 7i is isomorphic to 

nez nez 

where h(") := (e^"'^*) and lo = ^. 

According to Floquet theory of quantum mechanical systems driven by 
periodic perturbation [Ho,Yal,Ya2], the semi-group generated by K* is given 
by 

(e-^'^^7)(t) = U,{t,t- a) fit - a) , (58) 

where / G Til and a G R. Relation (j58|) can be seen by differentiating both 
sides with respect to a and setting a = (see [Ho]). (Alternatively, one may 
use the Trotter product formula, [RSI].) 
Note that if 

= , (59) 
for (j) = (j){t) G 7i and A G C, then (pit) satisfies 

[/^(t, 0)0(0) =e-*^V(^) • (60) 

Conversely, if 

^/^(r, 0)0(0) = e"^^>o, (61) 

then 

0(t)=e^"*?7,(t,O)0o (62) 
is an eigenfunction of K* with eigenvalue A. 



19 



We now study the spectrum of K* using complex spectral deformation 
techniques as developed in [HP,JP1,2,3].^ 

Let Uj be the unitary transformation generating translations in energy 
for the i^^ reservoir, z = 1, ■ ■ ■ , n. More precisely, for fi G i^^(R; B), 

MO)h{u) = f!{u) = f,{u + e). 

Moreover, let 

um :=r,(u,(e)) 

denote the second quantization of Uj(6'). 

Explicitly, Ui{6) = e^^^\ where Ai := idViidu^ is the second quantization 
of the generator of energy translations for the i*'*reservoir, z = 1, ■ ■ ■ ,n. We 
set 

U{Q) ■=\^ ®U^{Q)®...® Un{e) . (63) 

Let 

K;{e) := u{e)K;u{-e) (64) 
= -tdt + L;it,e) , (65) 

where L*(t, 6) is given by 

L;(t, 6) := U{e)Ll{t)U{-e) =U + Ne + gV^^^it, 9) , (66) 
Lo = L^ + J2i = dT{ui), z = 1, ■ ■ ■ , n, and 

i 

- zl^ ® a„ ® (-l)'^'(ai(e^'(">-^')/V£l(t)) 

- zl^ ® ® (-l)^'a*(e-ft(--'^»)/V£l(t))} • 

It follows from assuption (A2) that, for 6 G /(5), Vg°^{t,6) is a bounded 
operator. Hence K*{t,6) is well-defined and closed on the domain T> : = 
©(z^i) nV{N)n P(L^i) n ■ ■ ■ n P(L^"). when the coupling ^ = 0, the pure 
point spectrum of Lq is crpp(Lo) = {—2uo, 0, 2a;o}, with double degeneracy at 
0, and the continuous spectrum of Lq is crcont{l^o) = R- It follows that 

appiKo) = {E^''\g = 0) = + fcu; : j = 0, ■ ■ ■ , 3, A: G Z}, (67) 
^One may alternatively use the methods developed in [MMS1,2]. 
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where Eq i = 0, i?2 = —2uq, and = 2ujo, and (Jcont(-^o) = R- Let 

■= -idt + L^. 

Clearly, a{K^) = app{KQ). We have the following two easy lemmas. 

Lemma 4.1 

For 6 E C, the following holds. 

(i) For any tp G V, one has 

\\Ko{e)^^ = \\Ko{Ree)ijf + |/m^n|Ar^||2 . (68) 

(ii) If ImO 7^ 0, then Kq{9) is a normal operator satisfying 

Koiey = Ko{9) , (69) 

and V{Ko{e)) = V. 

(Hi) The spectrum of Kq{6) is 

(Tcont{Ko{0)) = {ne + s:ne N\{0} and s G R}, (70) 
app{Koie)) = {kuj + : J = 0, ■ ■ ■ , 3, A: G Z}, (71) 

where Eq i = 0, E2 = —2ujq and E^ = 2luq, (the eigenvalues of iP), and 

T 

Proof. The first claim follows directly by looking at the sector where = nl, 
since Kq{6) restricted to this sector is reduced to 

Kl''\d) = K^ + si + --- + sn + n9 , (72) 

where Si, ■ ■ ■ , s„ are interpreted as one-particle multiplication operators. For 
ImO 7^ 0, it also follows from (f?^ that 

V = {ij = : G P(<"V)) anrf ^ || J^'V)^^"^!!' < 00}, 

n 

and hence Kq{9) is a closed normal operator on T). Claims (ii) and (iii) follow 
from the corresponding statements on K^\9). □ 

Lemma 4.2 

Suppose (A1)-(A3) hold, and assume that {g,6) G C x I~{5). Then the 
following holds. 
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(i) V{K;{9)) = V and {K;{9)y = K^{9). 

(a) The map {g, 9) Kg{9) from C x I~{6) to the set of closed operators 

on 7i is an analytic family (of type A) in each variable separately; (see 
[Kal], chapter V, section 3.2). 

(Hi) For finite G R and Imz large enough, 

s - lim - zY' = {K;{Re9) - z)"' . (73) 



Proof. The first claim (i) follows from the fact that gV^°^{t,9) is bounded 
for 9 G It also follows from assumption (A2) that {g,9) K*{9) is 

analytic in ^ G Analyticity in g is obvious from ()65p. We still need to 

prove claim {Hi). Without loss of generality, assume that Re9 = 0. It follows 
from assumption (A2) that the resolvent formula 

iK;{9) - z)-' = {iq{9) - z)-\l+gV'"\;9){K*{9) - zyY' , (74) 

holds for small g, as long as z belongs to the half-plane {zGC:0<c< 
Imz}. Since {Kq{9) — z)~^ is uniformly bounded as Im9 f for 5' G R finite 
and Imz large enough, and V^°^{9) is bounded and analytic in 9, claim (iii) 
follows from the Neumann series expansion of the resolvent of K*{9). □ 

Next, we apply degenerate perturbation theory, as developed in [HP], to 
compute the spectrum of K*{9). Using contour integration, one may define 
the projection onto the perturbed eigenstates of K*{9), for Im9 G I~{S). Let 

P,m{0)--=j^ ^(^-A7(^))-\ (75) 

where 7^ is a contour that encloses the eigenvalues E^^\g = 0), j = 0, ■ ■ ■ , 3, A; G 
Z, at a distance > 0, such that, for sufficiently small \g\ (to be specified 
below) the contour also encloses Ej''\g). 

Moreover, let := -Po,(fc)-Pg,{A:)(^)-Po,(fc)-^° We show in Theorem 4.3 that 
the isomorphism 

4(fc)(g) :=T;(V)^Po,(fc)4(fc) (g) : Ran{Pg,^k)m ^ h« ^H^'^H'' (76) 

^°Note that, although Kq{9) is unbounded, it is a normal operator, and hence Po^(k) is 
well-defined by the spectral theorem; (see for example [Kal]). 
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has an inverse 

■= P9XkmPo,ik)f-ll)it) : hC^) ^n^c^n^-. Ran{P,,^^,^{e)). (77) 

We set 

M,,W := Po,wP,,w(^)i^;(^)P,,w(^)Po,w , (78) 
and define the quasi-Floquet Liouvillean by 

:= ^.,w(^)p.,w(^)i^;(^)p.,(.)(^)^,lfe)(^) = . (79) 

Let K = min{6, - ■ ■ ; ^l"}; where 6 appears in assumption (A2), section 
3, and Pi, - ■ ■ , Pn, are the inverse temperatures of the reservoirs T^i, ■ ■ ■ , Tin, 
respectively. For 6 G (see (I5ip ). we choose a parameter u such that 



-K < u <0 and - k < ImO < _!1iLM. (gg) 

We also choose a constant gi > such that 

giC <{k-\u\)/2, (81) 

where 

C := sup ||y*°*(t,^)|| < oo. (82) 

Theorem 4.3 

Suppose that assumptions (Al)-(A3) hold. Then, for gi > satisfying 
iSl]} . 9 G and v satisfying l\8(J\} . the following holds. 

(i) If \g\ < gi, the essential spectrum of the operator K*{9) is contained 
in the half-plane C\S(i/), where H(z/) := {2; G C : Imz > u}. More- 
over, the discrete spectrum of K*{9) is independent of 9 & U 
\g\ < ^gi, then the spectral projections Pg^{j.)[9), /c G Z, associated to the 
spectrum of K*{9) in the half-plane H(i^), are analytic in g and satisfy 
the estimate 

l|P,,wW-^o,(fc)|| <1. (83) 
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(a) If \g\ < y, then the quasi-Floquet Liouvillean ^g.{k) defined in l{79\) 
depends analytically on g, and has a Taylor expansion 

CO 

i=i 

where 

Kf^^ := kuj + ,k eZ. (85) 
The first non-trivial coefficient in ^84\ ) is 

=\i ^(^(^)("^(" " + (" ~ > 

where ^^k){z) = Po,(fc)^*°*(^)(^ - Ko{e))-'V'°\d)Po,(^k). 
Proof, (i) The resolvent formula 

me) - z)-' = {K*{e) - z)-\i + gv'"\;e){Km - zy^y^ , (86) 

holds for small g and z in the half-plane {zgC:0<c< Imz}. We extend 
the domain of validity of (j86p by refining the estimate on gV^°'^{t, 9){Kq{9) — 



Note that 



||^V^*°*(-,^)(i^*(^)-.)-^|| < \g\C\\{Ko[e) - z)-'\ 

< \9\C- ^ 



dtst{z,ri{K*{e))y 



where C is given by (jH^ and rj^KQ^O)) is the closure of the numerical range 



G{u,e) := {z EC : Imz > u; dist{z,'r]{KQ{9)) > e}. 



of Kq. Fix gi such that it satisfies (j81|) . and choose e such that e > > 0. 
Let 



Then 

sup ||^y*''*(.,^)(K*(e)-^)-i||<M. 

z&G(u,e) 9l 



24 



If \g\ < gi, the resolvent formula (|86p holds on G^u, e), and, for m > 1, 

m-~l f \g\ \m 

(87) 

It follows that 

U G{u,e)Cp{K;{9)), (88) 

''■^ 2 

where p{K*{6)) is the resolvent set of K*{9). Moreover, setting m = 1 in 
(jHTfl . it follows that, for \g\ < gi/2, 

\\Pg,{k){d) — -Po,(fc)|| < 1, 
and hence Pg,(k) is analytic in g. We still need to prove the independence of 

a^{K;{e)) oieei-{K). 

Fix {goyOo) e C X I~{6) such that \go\ < gi. The discrete eigenvalues 
of K*^{9) are analytic functions with at most algebraic singularities in the 
neighbourhood of 6*0, since K*^{6) is analytic in 6. Moreover, since K*^{6q) 
and K*^{6) are unitarily equivalent if {6 — 9q) G R, it follows that the pure 
point spectrum of K*^{6) is independent of 6. 

{ii) Analyticity of Tg,(fc) directly follows from (i) and the definition of Tg^i^^y 
Since — 1|| < 1 for \g\ < gi/2, ^l"^ is also analytic in g. Inserting the 

Neumann series for the resolvent of K*{6), gives 



with 



(90) 

Similarly, 

oo 

Mgm = K'' + Y.g'Ml^, (91) 
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with 



Mill = f ^Az-K^)-'Pomy'"'i-^m^-K^^^^^ 

(92) 

The odd terms in the above two expansions are zero due to the fact that 
Po,(fc) projects onto the = sector. The first non-trivial coefficient in the 
Taylor series of ^g(k) is 

^ = M^-l(fli^K- + K-f;i^) (93) 

= lf ^xu^){z - K'^r' + - K^r%{^)) , (94) 



with 

□ 



a(^) = Poik)V^''\;e){z - A-(e))-v;°*(-,e)Pow. 



We explicitly compute the discrete spectrum of K*{6) to second order 
in the coupling constant. Let ei_2 G TiP denote the vectors of spin up and 
down respectively. Then the states in 7i corresponding to the eigenvalues 
Ef\g = Q),3 = Q,--- ,3,/ceZ, are 



bl = e""^ ® ei ® ei 

bl = e*'"^ 62 ® 62 

bl = e*'"^ (g) 62 ® ei 

bl = e*'"^ ® ei ® 62 





■ ■ ■ ® fi^" 




■ ■ ■ ® li^" 











where ^2''^' is the vacuum state in JF^'(L^(R; B)). 

We apply perturbation theory to calculate E^l^\g). We know that 

^(') = ^/ |^m,(.)V^*°*(^^)(^-i^o(^))-^V'*°*(^^)Po,{.)(^-i^g))-^ 

+ (^- J^g))Po,(fe)V^*°*(-,^)(^-i^o(^))-'V*°*(-,^)Po,(fc)} . (95) 
For /^^^ as in (j^^ . we let its Fourier transform be 

W(n,^) := - r dte-'^'fp^^{u,t). (96) 
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Similarly, for as in (IHU)). we let 



fLJ^^^) ■■= - / dte—-'f*(u,t). (97) 



(k) 

Consider first the nondegenerate eigenvalue £"3 . Applying the Cauchy 
integration formula to and using the facts that 

limine = VV—, and lim Jm = 7iS(x), 

e\o X — le X t\o X — le 

where W denotes the Cauchy principal value, it follows that 

n 

/m«0|,S3'('Vfc» = -7r^^||/;3„M.,rn(2a;o-(/i;-m)u;,u;)|||, 

meZ 1=1 

where <^ ■, ■ ^ is the scalar product on Ti. Therefore, 

||/ft,/.„m(M,^)||| 



mGZ i=l 



{k — m)uj — u 



Yl ll/ft,M.,™,(2^o - {k - m)u, u)\\l + 0{g''). 



meZ i=l 

Similarly, 



Ei'\g) = kuj-2uJo-g'J2YVV du 



2ujq — {k — m)uj — u 



ing"^ 



mGZ i=l 

Next we use degenerate perturbation theory to calculate Ej''\g),j = 0, 1. 
Applying the Cauchy integration formula to and using the definitions of 
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and /*^, it follows that 



Re « 0r , ^r^T » = -Re « ^T'ct^i:' » 



0,1 vfc(2)^0,l ^ _ D„ ^ ^1,0 vfc(2)^0,l ^ 

||/ft,M.,m(M,w)||| 



^- — 1 ^ R. 



mGZ 2=1 



2u;o — {k — m)uj — u 



Im « 0^ , » = -Jm « 0^, tfUr » 



= — TT 

mGZ i=l 

Therefore, 

E^^}{g) = ku + g''ao,i + 0{g^) , (98) 
where ao,i are the eigenvalues of the 2x2 matrix 

mGZ i=l ^ ^ 

By construction, Kge^'^'^^Vt = kue^'^'^^n and t/(^)e*'='^*fi = e*^'^*fi, so {koo : k e 
Z} are also isolated eigenvalues of K*{6), 9 G This can be seen by defin- 

ing the spectral projections corresponding to the real isolated eigenvalues of 
Kg{6), using the resolvent, and taking the adjoint to define the corresponding 
spectral projections for the real isolated eigenvalues of K*{6)}^ 

The vector ip = is the eigenvector corresponding the eigenvalue of 

Hence, 

Ej,'\g) = ku;, (100) 

n 

E[''\g) = ku~ 2mg' E ll/ft-/^..-(2^o - (k - m)u)\\l + 0{g') . (101) 

mGZ 1=1 



^^Suppose a is an isolated and real eigenvalue of a closed operator A. Then the spectral 
projection corresponding to a is 

p=^j^ {z-A)-'dz, 

where is a contour enclosing a only. Since a is real and isolated, it is also an eigenvalue 
of A* with corresponding projection P* . (Using Cauchy's integration formula, one may 
readily verify that A*P* = P*A* = aP* and that (P*)^ = P*.) 
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Note that due to assumption (A3), ImEj < 0, for j = 1,2,3, while 
ItjiEU'^ = 0}^We have proven the following corollary to Theorem 4.3. 
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Figure 1. Spectrum of K*{9) 



Corollary 4.4 

Suppose assumptions (Al)-(A3) hold. Then, for 9 G and \g\ < 

gi/2, where gi satisfies \81\] . K*{9) has infinitely many simple eigenvalues, 
{fcti;}fcez, on the real axis, where = 

We will use the results of Theorem 4.3 and Corollary 4.4 to prove that 
the true state of the coupled system converges to a time-periodic state. 



5 Convergence to time-periodic states 

The following theorem claims that, under suitable assumptions, the true 
state of the system converges to a time-periodic state with the period r of 
the perturbation. 

Choose K = min((5, ■ ■ ■ ? ^) (as in section 4), define h**^** := _D(e'^v^^^"^^), 
and let 0*'^**''^ be the *-algebra generated by &*(/),/ G h*^^*, and by 1^. 
Note that 0*^^*-^ is norm-dense in 0^. We define a *-algebra 

C := 0^ ® 0^1.*^^* ® ■ ■ ■ (g) 0^-*^^* , (102) 

-"^^Alternatively, one can use the Feshbach map (see [BFS]) to compute the perturbation 
of the discrete spectrum of K*{d), which gives the same result. 



29 



which is dense in 0. 

We make the following additional assumption. 

(A4) The interaction Hamiltonian V{t) belongs to C, for t G R. 

Let 

D := 1^ ® 1^ ® e-V^^i+^ ® . . . ® e-'-v/^^ , (103) 

where A-ji. = dr{iduj,i = 1, ■ ■ ■ ,n, is the second quantization of the gener- 
ator of energy translations for the i^'^ reservoir. The operator D is positive 
such that DQ = Q, and RanD is dense in Ti. This opertor will be used as 
a regulator in order to apply the complex deformation technique. We have 
the following theorem. We will show in the following theorem that the time- 
periodic state to which the real state of the coupled system converges is re- 
lated to a zero-energy resonance state given by ^^3,(0) := (1 ®-D)-Pg,(o)(l Cg)f^) . 

Theorem 5.1 (Convergence to time-periodic states). 

Assume assumptions (A1)-(A4) hold. Assume further that a ^ C. Then 
there is a constant gi > satisfying 1181]) . such that, for \g\ < gi/2, the 
following holds 

hm (fi,«7+*(a)l]) = {ng,^o),D-'alia)n) , (104) 

where Clg.{o) corresponds to the zero- energy resonance of the adjoint of the 
Floquet Liouvillean, K*, and D is given by MOS^) . 

Proof. First note that by using a Dyson series expansion, it follows from 
assumption (A4) and the fact that a G C that a* (a) G C, and hence a* (a)r2 G 

The remainder of the proof relies on the result of Theorem 4.3, Corol- 
lary 4.4, and equation (j^Hj) (section 4). It follows from (j^Hj) and the time 
periodicity of f E H = L'^{[0,t]) that 

{e-'^>^l ^ n){0) = Ug{nT,0){l®n){0) = Ug{nT,0)n . (105) 

Let 1 ® =: n G TtC, and D -.= 10 0. 

Without loss of generality, we assume that uj = ^ 2ljo; (if ^ = 20^0, the 
state of the system typically oscillates between two resonance states until it 
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finally converges to a time-periodic state; see remark 1). Using the dynamics 
on C(0, fi) and ^TU^i . it follows that 

lim {n, a7+*(a)fi) = lim {Ug{nT, a* (a)fi) (106) 
= lim((e-*^>"n)(0),«*(a)fi) (107) 

Using the regulator D and complex spectral translation, 

lim a7+*(a)fi) = lim ((Df/(-^)e-'^^«(^)""?7(^);D n)(0), D'^a* (a)fi) 

(108) 



/oo 
du{u + ir]- i^*(^))-ie-*("+*'')"" 
-OO 



n){0),D-'al{a)n). (109) 
We split the above integration into two terms, 

/oo 
du{u + ir]- i^*(^))-ie-*("+^'')"^ 
-OO 

n)(0),D-^a*(a)fi) (110) 
= lim {(DU{-9) / ciz(z - K*(^))-^e-^"'^"n)(0), Z^-^a* (a)fi) + 

/oo 
dM(M - 2(/i - e) - ir*(^))-ie-^("-*('^-^))""n)(0), 
-OO 

^-^^^(a)^), (111) 

where ?7>0,0<e<yU, and 7 is the contour enclosing the point spectrum 
oiK*g{e) only. 

Using the results of Theorem 4.3 and Corollary 4.4, the first term con- 
verges to a time-periodic expression, 

lim {(DU{-9) I dz{z - A'*(^))"^e"^"""n)(0), D~^al{a)n) (112) 

n^oo I 

= J2 liin ((:Df/(-^)^-f.)(^)e-'^-w(^)"-4,w(^)^)(0),^-'«*(a)fi) (US) 
= J2{(DP,,(k)mm,D-'al{a)n) (114) 



fcez 
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Let 

:= DP,,(fc)(^)(e'''^* ® n) , (115) 
and denote by ^ ■, ■ ^ the scalar product on H. Then 

DPg,^k){e)D =< (e^^"* ® 1]), ■ > . (116) 

Therefore, 

DPg^^k)iO)^ =< e''^"* ® fi, 1 ® > (117) 

= ^^5,(0)4,0, (118) 

where 6kfi is the Kronecker delta, and hence 

J2{(DPg,ikmO),D-Wgia)n) = {in,,(^o)m,D-'alia)n) , (119) 

where ^^g,(o) is the zero-energy resonance of the Floquet Liouvillean. The 
second term in (jlllj) converges exponentially fast to zero since 

/oo 
dn(n-z(/i-e)-K;(^))-^e-*("-^(^-^))""n)(0),D-^a*(a)(]) = Oie-^''-'^^) , 
-oo 

(120) 

where < e' < e < yu; (see [Rud], chapter 19). □ 
Remarks. 

(1) When u = — = 2ujo, the system exhibits the phenomenon of reso- 
nance: The state of the system oscillates between two resonances until 
it finally converges to the time periodic state corresponding to ^g{o). 
This can be verified by a second order time-dependent perturbation 
theory calculation; (see also [Ya2]). 

(2) By a standard argument, the result of Theorem 5.1 can be extended to 
any initial state normal to u (see for example [JP3,MMS1,2]). 

(3) Note that 

((O,,(0))(0),D-^a*(a)fi) = ((17,,(,))(0),D-^(a)e^'=-*n) , 

where ^g.{k) is the state corresponding to the /ecu-energy resonance of 
the adjoint of the Floquet Liouvillean. In other words, all /cw-energy 
resonances belong to the same class of time-periodic states. 

In the next section, we discuss the positivity of entropy production per 
cycle and Carnot's formulation of the second law of thermodynamics. 
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6 Positivity of entropy production 



We consider a small system coupled to two fermionic reservoirs at the same 
chemical potential fi, yet at two different temperatures Pi and P2, with /3i < 
^2- Together with assumptions (Al)-(A4), we assume that the perturbation 
is differentiable in t, for t > 0. The first reservoir acts as a heat source, 
and the second reservoir as a heat sink. We want to show that, after the 
true state of the system has converged to a time-periodic state, the entropy 
production per cycle is strictly positive. We first prove that the time-periodic 
state, which the true state of the system converges to, is not normal to the 
initial state. 

We introduce the standard Floquet Liouvillean, 

kg := -%dt + l.g{t) , (121) 

acting on the extended Hilbert space i^^([0, r]) ® 7i, with periodic boundary 
conditions in t, where Lg(t) = Lq + gVif) — gJV{t)J, is the standard Liou- 
villean. We study the spectrum of the standard Floquet Liouvillean using 
complex spectral translations. Since the proof of the following proposition is 
very similar to the analysis in section 4, we only sketch the main steps of the 
proof. 

Proposition 6.1 

Suppose assumptions (A1)-(A3) (section 3) hold. Then there exists a 
positive constant g2, such that, for \g\ < g2, the spectrum of the standard 
Floquet Liouvillean Kg defined in M21]) . a{Kg), is absolutely continuous and 

a{kg) = aac{kg) = R. 



Sketch of proof . 

Let U{9) as in section 4. We define the complex deformed standard 
Floquet Liouvillean by 

kg{t,e) ■.= u{e)kgU{-e) (122) 

= -tdt + Lg{t,e) , (123) 

^•^In this section, we will use standard results about von Neumann algebras and Tomita- 
Takesaki modular theory. We refer the reader to [BR], chapters 2.4 and 2.5 for an expo- 
sition of these results. 
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where 
and 



v'^^t.e) = u{e){v{t) - jv{t)j)u{-e) 

n 
i=l 

It follows from assumption (A2) that V'^°*'{t,9) is bounded for 9 G Let 

C := sup (124) 

teR,6»G/((5) 

For ^ G I~{k), choose v such that > z/ > — k and — k < /m0 < — (k + |z/|)/2. 
Choose > such that 

g2C<{K + v)/L (125) 

Then, using an argument which is similar to the proof of Theorem 4.3, one 
can show that, for \g\ < g2, the essential spectrum of Kg{6) is contained in 
the half-plane {2 G C : Imz < z/}, and that it discrete spectrum 

app{Kg) = {Ef{g) : k e Z,j = 0, - ■ ■ ,3}, 
where (to second order in perturbation theory) 

3„M„m(M,^)||| 



meZ i=l 



2luq — {k — m)uj — u 



iTxg'^ 



E ll/ft,M.,™(2^o -{k- mp, 00)111 + 0(^?'), 

meZ i=l 

while 

E^%) = kuj + g'ao,i + Oig^) , (126) 
where ao,i are the eigenvalues of the 2x2 matrix 
2 



meZ i=l 



]_ _Q-l3i{2uio~{k-m)ui-fj.)/2 
g/3i(2t^o-(fc-m)t^-/^)/2 

(127) 
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Note that, to second order in the couphng g, the discrete spectrum is below 
the real axis. The claim of the theorem follows by noting that 



IrnOjO 



for small real g and large enough Imz. □ 




Figure 2. Spectrum of Kg (6) 

This result is sufficient to show that the time periodic state cu^^ defined 
in Q, section 2, is not normal to the initial state u. Let g^ := min{5fi/2, g2}, 
where gi satisfies flHTj) (section 4) and g2 satisfies 



Theorem 6.2 

Suppose assumptions (Al)-(A3) (section 3) and (A4) (section 5) hold. 
Then, for \g\ < g^, the time-periodic state cu^^,, for s G [0, r), does not belong 
to A/1,, ie, uj'^^g is not normal with respect to u. 

Proof. 

First note that under the assumptions of this theorem, the results of 
Proposition 6.1 hold. In particular. Kg has no real eigenvalues. 
By construction, 

o = ^+,. (128) 

Since we are assuming that cj^^ G AC, there exists a unique vector 
in the natural positive cone V := {aJaVL : a E M.} associated to the pair 
{M.,uj), such that 

a.+,(a) = (fi+„afi+,), (129) 



35 



hi aeM}^ Now, (gni), (II2HD, and (IT^ imply that 

(I7,(r, 0)fi+„ aF,(r, 0)1]+,) = afi+,), (130) 

for a G tVI. Using the facts that 

[j,Ug{t,t')]=0, (131) 

where J is the modular conjugation and Ug is the propagator generated by the 
standard Liouvillean, J^'^^s = ^g,s since eV, J* = J, and JAiJ = Ai', 
it follows that ()130|) also holds for a G Ai'. Furthermore, we know that^^ 

UgV C V, 

and that the linear span of V is dense in Ti. Let P^, be the orthogonal 
projection onto fi^^, then P^, G TVI V A4'. Moreover, since vr = vr^ ® tt''^^ 
• • • ® TT''^", in the Araki-Wyss representation, is a factor (of type II Ii, see 
[ArWy]), ie, M n M' = {Cl}}^ Suppose that 

I7,(r,0)fi+, = Cifi+, + C2vl>, 

where \Ef G Ran{l — P^^)' and ci,2 are complex numbers to be determined. 
Choosing a = (1 — P^,) in p30p gives 

|C2|' = 0. 

Together with the fact that 

(Z7,(r,o)fi+„F,(r,o)f]+,) = = i, 

it follows that |cip = 1. This implies that there exists A G R such that 

F,(r,0)l]+, = e-^-fi+,. (132) 

For each fixed s G [0, r], we define 

Ut):=e'''Ug{t,0)n+ (133) 



-"^^For a proof of this statement, see [BR], Theorems 2.5.31 and 2.3.19. 
i^see Corollary 2.5.32 in [BR]. 

i^Using the isomorphism between J^'^i {L^) (g) • • • (g) J^'^" (i^) and J^{L^^ © • • • ® L'^J, 
one can proceed to show that is a factor as in the case of the Araki-Wyss representation 
for a single reservoir of free fermions. 
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then ips is an eigenfunction of Kg with eigenvalue A. (This can be checked 
by looking at Kg(f)s-) However, this is in contradiction with the result of 
Proposition 6.1, and hence 0;+^ ^ AC- □ 

We have the following result regarding the strict positivity of entropy 
production per cycle. 

Theorem 6.3 (Positivity of entropy production) 

Suppose assumptions (A1)-(A4) hold. Then the entropy production per 
cycle, after the state of the system has converged to a time-periodic state, is 
strictly positive, ie, 

AEnt := [ dtiu^t{6^{gV{t))) > 0. (134) 
Jo 



Proof. It follows from assumptions (A1)-(A4) and Theorem 5.1, that 
sup I r dt{u^,^^,^iSUgVm-^oaliSU9Vm}\<oo. 

TGR+ Jo 

Together with the result of Theorem 6.2, this implies that u;^^ satisfies the 
assumptions of Proposition 2.1, and the entropy production per cycle after 
the state of the coupled system has converged to a time-periodic state, is 
strictly positive. □ 

Regarding the explicit computation of entropy production per cycle, fll34p . 
we would like to make the following remark. Since it follows from Theorem 
4.3 that i^g,(o) is analytic in g ioi g < gi/2, one can expand uj^^^ to any or- 
der in the coupling (7, and compute an explicit expression for AEnt given 
in ()134p up to this order in the coupling constant; (see also [MMS1,2] for a 
discussion of a perturbative approach for calculating entropy production in 
nonequilibrium steady-states). 
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Appendix A 



Glued Hilbert space representation 



We want to show that 



^(L2(R+; B)) ® ^(L2(R+; B)) = J^{L\R; B)) . 



Let VL be the vacuum state in the fermionic Fock space JF(L^(R+; S)). For 
fermionic creation/annihilation operators on JF(L^(R"'"; S)), 



define the creation/annihilation operators on JF(L^(R+; B)) ® JF(L^(R+; B)) 

as 



where ~ corresponds to complex conjugation. Note that hi and hr anti- 
commute. Let d and a* be the annihilation and creation operators on the 
fermionic Fock space JF(L^(R+; B) © L^(R"''; B)), such that they satisfy the 
usual CAR, and let VL be the vacuum state in {L"^ {R+ ; B) © L2(R+;i3)). 
An isomorphism between J^(L2(R+; B))®J^{L'^{Ii+- B)) and J^(L2(R+; i3) © 
L^(R+; B)) follows by the identification 



Now we claim that J^(L2(R+; i3)©L2(R+; B)) is isomorphic to J^(L2(R, rfw; S)). 
For 0, G R, consider the mapping 

J^,^ : L2(R+;i3) ©L2(R+;S) 3 (/, ^) ^ g L2(RrfM;i3) , 




hfif) :=6#(/)©l; 
6f(/) := (-1)^©6#(7), 



= «*((/, 0)), 
6f((7)-a#((0,<7)), 



such that 



h{u, a) : 



e*'''iy/m(|M|)5((|M|, (t) ,m < 



38 



This mapping is an isometry, since 

ll^llL2(R,dM;B) = Wif^ 9)\\l2{R+;B)<SL2(R+;B) 



dudam{u)\f{u, a)\ + / dudam{u)\g{u, a)\ 
n+;B Jb.+ ;B 

2 1 1 1 1 2 

L2(R+;B) + II5'IIl2(R+;S) • 

Moreover, the mapping j^^^p is an isomorphism, since, for given h G L'^(R,', B), 
there exists a mapping '■ h ^ (f, g) E L'^{IV"; B) © /.^(R"*"; B), such that 

0-) ^^^ ^(^, 0-), n > , 
^Jm{\u\) 
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